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Abstract 

In this paper "we present the time resolution measurements of the Lutetium- 
Yttrium Oxyorthosilicate (LYSO) calorimeter prototype for the Mu2e experi¬ 
ment. The measurements have been performed using the e~ beam of the Beam 
Test Facility (BTF) in Frascati, Italy in the energy range from 100 to 400 MeV. 
The calorimeter prototype consisted of twenty five 30x30x130 mm^, LYSO crys¬ 
tals read out by 10x10 mm^ Hamamatsu Avalanche Photodiodes (APDs). The 
energy dependence of the measured time resolution can be parametrized as 
at{E) = a/yjE/GeV 0&, with the stochastic and constant terms a = (51 ±1) ps 
and b = (10 ± 4) ps, respectively. This corresponds to the time resolution of 
(162 ± 4) ps at 100 MeV. 
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1. Introduction 


The Mu2e experiment at Fermilab [T] aims to search for Charged Lepton 
Flavor Violation (CLFV) in the neutrinoless, coherent conversion of a negative 
muon into an electron in the Coulomb field of an ^^Al nucleus. The /x —> 
e conversion results in monoenergetic electrons with an energy equal to the 
muon rest mass minus the corrections for the nuclear recoil and the binding 
energy of the muon. For ^^Al, the energy of the monoenergetic electron is Ef. 
= 104.97 MeV [5]. 

The experiment is designed to reach the single event sensitivity (SES) of 
2.4 X 10“^^ in three years of running [1]. This value represents an improve¬ 
ment of four orders of magnitude over the current best experimental limit 
i?^e(A-u) < 7x10“^^ @ 90% C.L. set by the SINDRUM II experiment [3]. 

The Standard Model predicted rate for this process is 0(10“®^) [3], therefore 
any signal observed by Mu2e would be a compelling evidence of new physics. 

2. The Mu2e electromagnetic calorimeter 

The Mu2e detector is designed to identify /x —)■ e conversion electrons and 
reduce the background to a negligible level. The detector is located inside a 
large superconducting solenoid with the magnetic field B=1 T and surrounded 
by the cosmic ray veto counters. 

A low mass straw tracker provides an accurate track momentum measure¬ 
ment necessary to separate the signal from the background. The calorimeter is 
located behind the tracker and complements it by providing: (i) powerful fi/e 
particle identification (PID), (ii) seeds for the pattern recognition in the tracker 
and (iii) an independent software trigger system. 

Efficient PID requires the calorimeter to have the timing resolution better 
than 500 ps, the energy resolution of 0(5%) is needed to provide an efficient 
trigger. The calorimeter should be able to operate in an environment where a 
radiation dose up to 120 Gy/year is delivered by protons, neutrons, and photons. 
It must also function in a I T axial magnetic field and a 10“^ Torr vacuum. 
Before the sudden increase of the lutetium price, the Mu2e calorimeter design 
included two disks of LYSO crystals read out by two large area avalanche pho¬ 
todiodes (APDs) per crystal [5]. The choice of LYSO as a scintillator provided 
high light yield, fast response, and radiation hardness [^{5]. 

In this paper, we report the results of tests performed with a LYSO-based 
calorimeter prototype, which include measurements of the timing resolution and 
evaluation of the front-end electronics (FEE) and readout system. 

3. The LYSO crystal calorimeter prototype 

The calorimeter prototype consisted of 25 LYSO crystals (30x30x130 mm^) 
from the Shanghai Institute of Ceramics, Chinese Academy of Sciences (SIC- 
CAS) [S], arranged in a 5x5 matrix (Fig. [^. Each crystal was identified by 
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Figure 1: On the left (right), is the front (back) view of the 5x5 LYSO calorime¬ 
ter prototype. The APDs that are attached to the back of each crystal are 
visible in both views. The Fig. on the right shows a brass Faraday cup that is 
placed around placed around each Amp-HV board. 


two indices (i,j), corresponding to its row and column positions, with the crys¬ 
tal (0,0) being the bottom-left crystal in the matrix viewed from the back. 
The dimensions corresponded to ^ 11.2 radiation lengths (Xq) in depth and a 
transverse size of ~ 3.6 Moliere radius (R^f). The crystals were individually 
characterized with a ^^Na source and a spectrophotometer, 15 crystals at INFN 
Laboratori Nazionali di Frascati (LNF) and 10 at the California Institute of 
Technology (Caltech). All crystals demonstrated high light yield and transmit¬ 
tance. The measured longitudinal response uniformity was below 5%. Each 
crystal was wrapped in a 60 pm thick layer of a super-reflective 3M Enhanced 
Specular Reflector (ESR) film [TU] and read out by a Hamamatsu S8664-1010 
APD m- The APDs were optically connected to the crystals with Saint-Gobain 
BC-630 optical grease [T^ . 

The Front-End Electronics (FEE) consisted of a multi-layer, double-sided, 
discrete board (Amp-HV) directly connected to the photosensor pins (see Fig.[^. 
The board provided both the amplification stage and the local regulation of the 
photosensor bias voltage, thus reducing the noise loop-area. The amplification 
layer was a double stage transimpedance preamplifier with a total gain of 15 kD, 
which maintained an equivalent noise charge (ENC) level of about 10^ electrons 
with no input capacitance source. The linear regulation layer allowed precise 
voltage regulation and long-term stability of better than 100 ppm. Each group 
of 16 Amp-HV chips was controlled by an external ARM controller distributing 
the voltages. Two ARM controllers were used for the prototype. The high 
voltage (530 V) was produced by a primary generator that used the low-noise 
switching technology and resided on the ARM controller board. The output 
voltage was regulated by a DAC and read out with a 16-bit ADC. 

APDs were calibrated using a green (530 nm) 0.6 pJ/pulse laser. The laser 
light was distributed through 250 pm core diameter fused silica optical fibers. 


3 




Figure 2: A picture of an Amp-HV prototype. Left: amplification. Right: HV 
side. 


The fibers were inserted into a custom connector, polished and positioned di¬ 
rectly in the APD holders. The 10 ns laser pulse was synchronized with an 
external trigger. 

4. Experimental setup at BTF 

The calorimeter prototype was tested at the Beam Test Facility [13] (BTF) 
of the INFN LNF in December 2014 and again in April 2015. The BTF is a 
part of the DA<i>NE (Double Annular Factory for Nice Experiments) accelerator 
complex equipped for testing particle detectors. The complex includes a Linac 
which sends the beam pulses to the BTF area at a 50 Hz frequency. Each 
Linac pulse has a ~10 ns duration and is divided into 180-200 ps long bunches. 
The trigger signal provided by the Linac has a time resolution on the 0(10 
ns) time resolution, making it necessary to use a different source of the trigger 
timing. Two 6x10x50 mm^ finger-shaped beam scintillation counters located 
upstream of the calorimeter prototype have been used for this purpose. Due 
to interference with other detectors and logistic in the area, we were unable to 
place the counters closer than 60 cm to the calorimeter surface. The scintillation 
counters were read out by 3x3 mm^ SensL [14] silicon photomultipliers (SiPMs) 
(Fig.i- 

In order to select cosmic rays, two plastic scintillation counters, 50x50x200 
mm^ in size, were positioned above and below the calorimeter prototype. Each 
of those cosmic counters has been read out by two photomultipliers (PMTs). 

Initial channel-to-channel calibration was performed with cosmic ray min¬ 
imum ionizing particles (MIPs) and later updated using more accurate beam 
calibration. 

The data taking configuration used an OR of three different triggers: (1) a 
beam trigger (BT) formed by the AND of signals from the finger scintillation 
counters. The rate of this trigger varied from run to run, from a few Hz up to 
20 Hz; (2) a cosmic ray trigger (CRT). The rate of this trigger was at a level of 
2 Hz; (3) a laser trigger (LT) generated by a timer at a typical rate of 0.1 Hz. 
The CRT and LT triggers were also used to monitor the calorimeter response 
during the data taking. 
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Figure 3: Left: picture of the BTF hall with the experimental setup used for the 
calorimeter prototype time resolution measurement. Right: schematic diagram 
of the experimental setup seen from the top. 


Data were acquired by a VME-based DAQ system from CAEN m and read 
out by 4 CAEN VI720 waveform digitizers at a sampling rate of 250 Msps with 
12-bit resolution over the 0-2 V dynamic range. 

For each calorimeter channel and each trigger, the signal charge was deter¬ 
mined by numerical integration of the signal waveform in a gate of [-50, 450] ns 
around the time sample corresponding to the pulse maximum. The charge base¬ 
line was estimated by integrating the waveform in a gate [-750, -250] ns before 
the pulse maximum. The baseline was then subtracted from the signal on event 
by event basis. 

5. Equalization and calibration 

In order to measure the time resolution of the calorimeter prototype in the 
energy range from 100 to 400 MeV, two different APD gain settings were used. 
A gain of G = 75 was used for runs at 100, 150 and 200 MeV. For signals to stay 
within the dynamic range of the digitizer a gain of G = 25 was used at higher 
beam energies. The laser trigger was used to adjust gains in individual channels. 
Waveforms from each channel were monitored with a software scope and the 
pulse heights were equalized by adjusting the individual APD HV settings. 

The crystal-to-crystal response was determined by directing 450 MeV elec¬ 
trons onto the center of each crystal. The global charge-to-energy conversion 
scale was set by comparing the total energy reconstructed in the calorimeter, 
Erec, to the expected energy deposition estimated by the Monte Carlo (MC) 
simulation based on GEANT4 m- 

6. Event selection 

The BTF beam intensity was tuned to provide the mean number of electrons 
per bunch A ~ 0.8 at 100 MeV. At this intensity a non-negligible fraction of 


5 











bunches contains two or more electrons. In order to select only single-electron 
events, a cut on the total reconstructed energy in the calorimeter prototype, 
Erec^ and the signal charge in each of the scintillation counters, Qdep, has been 
applied. The cut on the energy deposited in the calorimeter prototype has been 
set to Free < 1.3 • Ebearm as shown in Fig. I^left). An example of a Qdep cut is 
shown in Fig. gright). 




Figure 4: Left: reconstructed energy in the calorimeter prototype for an e~ 
beam of 100 MeV. Right: charge response for one of the scintillator counters 
for an e~ beam of 100 MeV before and after the cut on the energy deposited 
in the calorimeter. The single-electron selection cut has been set to A^ec < 
1.3 • Ebeam for the reconstructed energy in the calorimeter prototype and to 
80 pC < Qdep < 160 pC for the charge deposition in the counter, as shown by 
the dashed lines. 


6.1. Centroid cut 

Because of ~60 cm distance between the calorimeter prototype and the scin¬ 
tillation counters, the effect of the electron multiple Coulomb scattering down¬ 
stream of the counters was not negligible, especially at low energies. In Fig. 
the correlation between the x (horizontal) and y (vertical) coordinates recon¬ 
structed in the calorimeter for Ebeam = 100 MeV is shown. The x {y) coordinate 
is given by the logarithmic energy-weighted average of the crystal positions 

25 

Y, Xi [wo - \og{E^/Etot)\ 

-- ( 1 ) 

Y ko - ^og(EjEtot)\ 

i=l 

where Xi can be 0 cm, ±3 cm, ±6 cm (centers of the crystals with respect to 
the beam line), Ei is the energy deposition in the i-th crystal, Etot is the total 
energy deposited in the calorimeter and Wg is a custom parameter set to 9. To 
suppress effects of multiple scattering we require 

r = ■\/x^ + y^ < 0.5 cm. (2) 
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Figure 5: Spatial distribution of hits on the front face of the calorimeter proto¬ 
type at Eteam = 100 MeV. The dashed lines show 9 of 25 crystals. 
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As the multiple scattering scales as 1/p, at higher beam energies its effect be¬ 
comes smaller. 


7. Determination of the time resolution 

To extract the calorimeter timing information, the waveform signals asso¬ 
ciated with e~ events have been fitted with a Landau function in the range 
[tmax — 30 ns, tmax + 70 ns], where tmax is the time of the waveform peak, as 
measured by the digitizer (Fig.j^. Signals from the trigger counters have been 
fitted with a log-normal distribution [17] in the range [tmax —15 ns, tmax -f 10 ns]. 




Figure 6: Examples of the calorimeter signals. Left: waveform for a e~ beam 
event fitted with the Landau distribution. Right: waveform for a laser event 
fitted with the log-normal distribution. 


The time corresponding to the maximum of the fitted function has been used 
as the signal time. Signals above 10 mV were used in the data analysis. 

7.1. Calorimeter time resolution with an external .start 

The time resolution was determined from the width of At = tgignai — tstart 
distributions, using either the central crystal alone or the entire calorimeter 
prototype. In the latter case, tsignaZ is an energy-weighted sum of times recon¬ 
structed in different channels: 


25 

z; uEi 


t 


signal 




25 
2 = 1 


( 3 ) 


where C is the peak time of the z-th crystal, is the energy deposition in the 
z-th crystal, and tstart = {tf^ Etf.^)/2 is the average of the beam counter times. 

After correcting for delays due to the cable length differences between chan¬ 
nels (T°s), a residual time-walk effect, a dependence of U on the deposited 
energy, remains. In Fig. examples of this dependence for the central crystal 











for Ebeam = 100 MeV and Eteam = 400 MeV runs are shown. Due to the dif¬ 
ferent gains (25 vs 75), the energy depositions at Eteam = 400 MeV run have 
been scaled down by a factor of 3. We parameterize the dependence with a 
a + b/E + c/\/~E function and correct the individual reconstructed times in Eq. 
<§ as follows: 

, bi Ci , . 



Figure 7: Central crystal timing (t 22 ) as a function of the deposited energy 
in the central crystal for Ebeam = 100 MeV and Ebeam = 400 MeV. Energy 
deposition for Ebeam = 400 MeV has been corrected by a factor of 3 due to the 
different gains (75/25). 


Fig. 1^ shows distributions of the corrected tgignai at 100, 200, 300, and 400 
MeV with the gaussian fits superimposed. Gaussians describe well the central 
part of the distributions, a small contribution of 2-electron events leads to non- 
gaussian tails. By observing the relative shifts of the means of the distributions 
with respect to zero, the precision of the calibration procedure was estimated 
to be better than 20 ps. 

The dependence of the time resolution on E^ep is shown in Fig. and is well 
described by the function 
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(c) 300 MeV 


(d) 400 MeV 


Figure 8: Corrected distributions of tmatrix for 100, 200, 300, and 400 MeV. 
The long tail present on both the high and low sides is related to a small 
contamination of two-particle events that are still present in the selected sample. 
However, the high-side tail in the time distribution at 200 MeV (top-right) is 
related to the pulse height saturation observed at that energy. 
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where a is the stochastic term and b is the constant term, which is mainly due 
to the trigger jitter. From the fit, a = (51 ± 3) ps and b = (157 ± 7) ps. 

The trigger jitter can be estimated directly by htting the distribution of 
with a gaussian, as shown in Figure 10 The width of the distribution 
returned by the fit is cr = (287±5) ps. Assuming the resolution of both counters 
is the same 

cr {tstart) = -th) = 144 ± 3ps, (6) 

consistent with the constant term from the energy dependence fit ©■ 



Figure 9: Time resolution, as defined by the central crystal and the energy- 
weighted sum over all crystals, as a function of the deposited energy Edep- The 
dashed line represents the stochastic term only. 


1.2. Calorimeter-based time resolution 

The calorimeter time resolution can also be determined by measuring the 
time difference between the signals in the two neighboring crystals. This tech¬ 
nique does not require an external time reference and is widely used in HEP |18] . 
This method has been applied to the data collected at 100 and 200 MeV with 
the beam offset horizontally by 0.6, 1.0 and 1.5 cm with respect to the prototype 
center. Events with the reconstructed energies in the neighboring crystals E 12 
and E 22 satisfying the requirement: 

0.8 < E 12 /E 22 < 1-2, (7) 
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N. Entries / 40 ps 



Figure 10: tf^ ~tf 2 distribution for a 100 MeV e beam. The half width of the 
gaussian gives an estimate of the trigger jitter. 
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have been selected for analysis. Time-walk corrections in individual channels 
have been determined as described in section ITTI 



Energy [MeV] 


Figure 11: Time-walk corrections for the run at 100 MeV. ^22 is the crystal tim¬ 
ing measurements for the central crystal before the correction. The fit function 


is described in section 7.1 


Fig. [TT] sh ows an example of the time-walk correction fit for Ebeam = 100 
MeV. Fig. |12| shows the ti 2 —122 distributions for beam energies of 100 and 200 
MeV. 

7.3. Time resolution with MIPs and laser 

Measurements using minimum ionizing particles (MIPs) allow the time res¬ 
olution at energies significantly below 100 MeV to be determined. For that, 
we used data collected with the cosmic trigger described in section]^ Selected 
events were required to have the energy deposition above 5 MeV in a column of 
5 crystals, surrounded by two columns with energy depositions below 5 MeV in 
each. 

The time resolution has then been measured for a single crystal, using as 
tstart the time of another crystal of the same column. The time resolution 
averaged over multiple tested pairs of crystals was (312±9) ps for a mean energy 
deposition of (23.1 ± 0.3) MeV, estimated using the GEANT4-based simulation. 

The same procedure has been applied for the energy-weighted time average 
of two crystals, using as tstart the time average of the two other crystals in the 
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Figure 12: Time residual distributions between signals in two adjacent crystals 
for Ebeam = 100 MeV (left) and Ebeam = 200 MeV (right). 


same column. The time resolution corresponding to the mean energy deposition 
of (46.2 ± 0.6) MeV was (262 ± 7) ps. 

In order to estimate the contribution to the timing uncertainty from the 
photosensor, FEE, and digitizer, the data collected with the laser trigger at the 
APD gain G = 75 have been used. 

Ten out of 25 crystals have been chosen to provide a reference time 

10 

t start = -■ (8) 

1=1 


The rest fifteen crystals have been split into 5 groups of 1, 2, 3, 4, and 5 
crystals correspondingly. For each group, an energy-weighted mean time has 
been calculated as 


tr), — 


i=l 

n 


(9) 


and the widths of the distributions in — tgtart have been plotted versus the 
mean total energy corresponding to the integrated charge in the crystals in¬ 
cluded in the group, as shown in Fig. 13 A ht with the function: 


(Jt{E) = 0 5, 

v/A/GeV 


( 10 ) 


returns a = (9 ± 1) ps and b = (48 ± 4) ps. 


8. Conclusion 

The calorimeter time resolutions measured using different techniques are 
shown together in Fig. Measurements with the external start are corrected 
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Figure 13: Time resolution for laser signals as a function of the equivalent 
deposited energy. 
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for the timing jitter of the start signal. Measurements using the neighboring 
crystals method are converted to the single channel timing resolution assuming 
the resolution of both channels is the same. 

The energy dependence of the timing resolution is parameterized as 



( 11 ) 


with the stochastic term a = (51 ± 1) ps and the constant term b = (10 ± 4) 
ps determined from the fit. 

The time resolution of the LYSO calorimeter prototype at 100 MeV, at = 
(162 ± 4) ps, amply satisfies the Mu2e calorimeter requirement at < 500 ps [T]. 

The time resolution also can be presented as the quadratic sum of three 
terms 


( 12 ) 


0 O'fee 0 O: 


Ot = 


VA^p.e./MeV 


where the first term is due to the photo statistics and the emission time of 
the scintillator the second term corresponds to the timing jitter due to the 
photosensor and electronics, and the last term, ax, accounts for the shower 
length fluctuations, reconstruction and calibration-related terms. 

For Ts(LYSO) = 40 ns [19] and the number of photoelectrons estimated at 
Np.o./MeV = 4100, the first term at 100 MeV contributes around 63 ps. Ac¬ 
cording to Eq. pr] ), at this energy, a fee = 38 ps. The remaining contribution, 
ax ~ 144 ps, is dominated by the waveform reconstruction technique used. 

Due to the limited energy coverage of our data is difficult to estrapolate 
this results to higher energy range (> GeV) where the stochastic part of this 
dependence will be practically negligible. In order to reach timing close to the 
constant term i.e. 0(10) ps, there will be need to improve both the front end 
electronics, that should demonstrate to have a rise time independent from the 
pulse height, and the fitting technique of the pulse shape. 
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